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I. INTRODUCTION 



In a recent paper we presented a simple, formally covariant scalar-diquark-quark 
model of the nucleon [|l|. Combining the Salpeter approach with the Mandelstam 
formalism we calculated the electromagnetic form factors of the nucleon. Despite 
the simplicity of the model, we obtained a very good description of the measured 
proton electric form factor for momentum transfers up to —3 GeV 2 . However, the 
electric form factor of the neutron came out a factor of more than three too large, 
and the calculated magnetic form factors failed to describe the data quantitatively. 
Since scalar diquarks do not contribute to the magnetic current density, the magnetic 
form factors result only from the coupling to the quark. This is the main motivation 
for introducing also spin 1 diquarks. Of course, in the non-relativistic quark model, 
two quarks may couple to spin or to spin 1, thus forming a mixed-antisymmetric 
or mixed-symmetric spin state. One main idea of the present model is to identify 
the mixed-antisymmetric part of the spin function with the scalar diquark channel 
and the mixed-symmetric one with the axial- vector diquark channel. An axial- vector 
particle has positive parity, and we consistently neglect any relative angular momen- 
tum between the diquark and the quark. There are indeed theoretical arguments 
in favour of a dominance of scalar diquarks @,fj, and also of eventual axial- vector 
diquarks [f|||. Also, the amount of experimental hints of diquarks in the nucleon 
is not negligible, see [[]] for a review on that subject. On the other hand, there are 
theoretical works which clearly deny such correlations, see e.g. J7|. However, we do 
not want to discuss against or in favour of diquark correlations. Obviously, this is 
only possible in a full three particle calculation. Rather, we develop a simple but for- 
mally covariant model of the nucleon which is able to reproduce various experimental 
data such as electromagnetic form factors, mean square charge radii and magnetic 
moments. 

As in |]J we start from the Bethe-Salpeter equation. Following Salpeter || we 
assume an instantaneous interaction and obtain a Salpeter-type equation. As the 
only interaction we adopt an (instantaneous) quark exchange between the quark and 
the diquark. This interaction has been previously used in various studies [p|-[l^j5[1 
and seems most natural for a quark-diquark model. Involving only scalar and axial- 
vector diquark channels, we deduce a pair of coupled integral equations similar to 
12fl . Within a basis of positive parity amplitudes with spin ~, we obtain a bound 



state solution below the threshold for this Schrodinger-type equation. Then, electro- 
magnetic transition currents and form factors are calculated using the Mandelstam 
formalism [[TJJ. For details of the calculation see |IJ. 

This paper is structured as follows. In Sec. [H] we summarize the basic equations 
of the model. The pair of coupled Salpeter equations is presented. The hermiticity 
condition for the interaction is shortly outlined. Sec. |HI| presents the idea of how 
the proton and neutron currents are calculated. The results of some current matrix 
elements and of the form factors are given in Sec. |V]. We compare with experimental 
data. Finally, a summary is given in Sec. |V|. 
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II. THE MODEL 



A. The diquark-quark Salpeter equation 

A relativistic bound state of a scalar or (axial-)vector particle and a quark with 
total four momentum P is described by a Bethe-Salpeter amplitude 

Xp(xi,x 2 )( fl ) a = (0 T 4> M (x 1 )i! a {x 2 ) P) ■ (1) 



In momentum space, xp fulfills the following Bethe-Salpeter equation [1TB 



Xp(p) W = A f (Pi)$SZ (Pa) / t^s HK(P )P ,p')XpO/)) M (2) 



with Af (pi)(^ and 5 , ^(p 2 ) the Feynman propagators for a scalar (vector) particle 
and a quark, respectively, and — %K the irreducible interaction kernel. The fundamen- 
tal equation of our model has been derived in 0. It is a Schrodinger-type equation 
of the following form (all indices suppressed): 

(HV)(p) = M V(p) 

= (r+w)*(p). (3) 

^(p) is the Bethe-Salpeter amplitude in the bound state's rest frame integrated over 



pO. 



! • (4) 

P=(M,0) 



W(p,p') = K(P, p,p')7° is the interaction kernel in the rest frame which is supposed 
to be independent of p° and p'°. H^ip) is the Dirac Hamiltonian of the quark, and 
= P 2 + mf the squared on-shell energies of the diquark and quark, respectively 
The masses are constituent masses of the order of a few hundred MeV, see Sec. 
[TV] . M is the mass of the bound state. For the normalization we obtained in |TJ : 

(* = 2M (5) 

with the scalar product defined as 

|¥> = / ^(2^) *%) *(p) . (6) 

To obtain real eigenvalues M of the pseudo-Hamiltonian 7i in Eq. (|3|) the interaction 
kernel has to be chosen such that 7i is hermitian with respect to the scalar product 
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B. The interaction kernel 



Scalar and axial-vector diquarks (called v-diquarks in the following) couple to two 
quarks according to the Lagrangians 



eg*-™** = - 9v v c r r a $ , (7) 

respectively. Here, g s and g v are dimensionless coupling constants, and we will choose 
9s — 9v = 9, see Sec. |V[ To obtain a stable solution of the Salpeter equation (|3|) 
(or (|), see below), we have to introduce a form factor of the diquarks, see It is 
assumed to be the same for the scalar and v-diquark and is chosen to be a Gaussian 
function. Again we use as the only interaction between the diquark and quark a 
quark exchange and perform the instantaneous approximation, i.e. neglect the p° 
dependence. So the interaction kernel is of the form 



W(p,p') = -g 2 -^(-j(p + p')+m q )j° 



q 



with uj q = \J(p + p') 2 + m 2 the energy of the exchanged quark. We redefine the 
above g by absorbing the colour factor of 2 (For a given colour of the diquark there 
are two possible colours of the exchanged quark). We then obtain from Eq. @ the 
following coupled integral equations for the scalar and v-diquark components of the 
nucleon, namely and ( see Fig. P (without flavour dependence): 



U)i + Co>2 
^2 

1 r d 3 p' 2 1 



M ¥°\p) = ^l^lH 2 {p)¥°\p) 

L0 2 



2wi J (2vr) 3V 9s 'uo, 



-l(p + p')+m q ) 1 Q ^(p') 



^(-9sg v )^lH-i(P + P)+rn q h°^ ] (p') 



1 

+ (2tt) 3V aaav 'u\ 



M^{p) = ^^H 2 {p)¥^{p) 
^ uj 2 



~ / 7|^(+^)^TV(-7(p + ^)+m 3 ) 7 5 7^7 *L 11 (p , ) 
1 r d 3 p' . . 1 



2u x J (2tt) 3V co 



(+gsg v )-J-i(p + p')+m q hHi°y m (p) • (9) 



The above interaction W (compare with Eq. @) is indeed hermitian with respect 
to the scalar product @, since the following conditions are fulfilled: 

www&py = w®w(p',p) 

wN^MfY = www\(p',p) (io) 
w^°^(p,p'Y = w^{f,p) , 
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where e.g. 




The sign of PI/' " ] is chosen in [|l| to obtain a bound state of only a scalar diquark 
and a quark. With respect to a further extension of our model (see Sec. [V]) we choose 
the sign of WWW to be positive. It can be shown that this sign gives rise to a bound 
A(3/2) composed of a v-diquark and a quark. 



C. Solving the diquark-quark equation 

The procedure solving Eq. (Q) is analogous to |1] and involves the Ritz variational 
principle. As a basis of the nucleon with spin s) in spin and Lorentz space we 
choose 

ef 1 = \el, 0, 0) , e; [ J = *(0, e* t 0) , = '(0, 0, e*) , 



with e s (p) = 











(s) 




) = 


V W+M12 ^ s / 



(12) 



where M 12 = mim2 and W 2 = p 2 + M, 2 ,. Note that we combined the three vector 
components of the v-diquark to one reduced function only carrying total spin 1. 
The corresponding basis vector is denoted e l }y- The contribution from the zero- 
component of the v-diquark is expanded in terms of e* ) and that from the scalar 
diquark channel in terms of e* . However, the basis states of Eq. ( |T2"D are not those 
involved when calculating the matrix elements of the potential in Eq. @. There, 
the 7 matrices are coupled with the quark spin. Thus, for the quark coupling to the 
zero-component of the v-diquark (via ~ 7 5 7°), also the following basis states with 
negative parity are needed (see App. |A|, |B]): 

e 3 s (p)= ( °J , et(p)= (w^X'^J (13) 

Then, 





7 5 7° e? s (p) 


~ e \ip) 




7 5 7°eKp) 




7 5 


7 [1] <g> e 1 (p) 


1 

2 

s 




7 5 


7W <g> e 2 (p) 


1 

2 
S 


~ e s(p) 



(14) 



are basis states of positive parity, see also App. 
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III. CURRENT MATRIX ELEMENTS AND FORM FACTORS 



The current matrix elements are calculated in the Mandelstam formalism [13 
The corresponding diagrams are shown in Fig. |2|. For details of the calculation 
see M. The Bethe-Salpeter amplitudes of course depend on the total and relative 
momenta P and p. A correct boost prescription (App. |B|) for the outgoing \ is 
crucial for a relativistic treatment. Since the implementation of this formalism in 
the scalar diquark-quark model is elaborated in [1]] we only focus our attention on 
aspects concerning the v-diquark. 
The propagator of the v-diquark is chosen to be 

9nu 



i.e. we neglect the term 



Af(p) 



p2 _ m 2 _|_ i e 



(15) 



There are several technical reasons for this: Firstly, 



it is not possible to state a coordinate independent normalization condition with the 
full propagator. Secondly, in the Salpeter approach, the contour integral over the full 
propagators in the complex plane does not vanish at infinity for the components with 
fi,v = 0. Finally, the full propagator has no inverse, which, however, is mandatory 
in the Salpeter approach, e.g. to derive the normalization condition. 
The coupling of a photon with index fi to a massive vector particle is given by p~T|JT8[1 
(see Fig. ||): 

IV; ba = ~(P + P% 9ba + ((1 + K-)Pb ~ (« + £)p'J Qua 

+((i + «K-(« + 0p«)&6, ( 16 ) 

with k the anomalous contribution to the magnetic moment of the v-diquark (k = 1 
for a pointlike spin-1 particle). £ is a gauge parameter introduced by Lee and Yang. 
According to our choice of the propagator ( |l5f) we put £ = 1. We will see that this is 
indeed necessary to obtain a conserved current and the correct normalization of the 
Salpeter amplitudes. 

Rewriting the coupling ( ]T5| ) in spherical components with indices a and b and taking 
into account the corresponding Clebsch-Gordan coefficients, we obtain the coupling 
matrices in the space of (e' ', e ey ), i.e. including both scalar and v-diquark cou- 
plings: 



f(pi+p'i)° 




V 



V 



o o \ 

-(Pi+Pi)° + 

l(l + „)g {p 1+ p[) 







-(pi + p\f 













(pi + p'l) 3 
-^(i + «)g° 





/ 

\ 

^(1 + K)g° 

-(Pi +P'i) 3 



A/ 1 



(17a) 



(17b) 



(17c) 
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where T + = ^(r 1 + iT 2 ). (q°,q) is the photon momentum and q — \q\ — — g 3 . 
Transitions between scalar and v-diquarks are neglected since they would change the 
flavour symmetry in our model [pT3| , p!8| . From Eq. (|17a| ) we see that, in accordance 
with the choice of the propagator, the zero component of the v-diquark contributes 
negatively to the norm. As in we see that the current is conserved exactly also in 
the v-diquark channel. 

In the non-relativistic quark-model the spin-flavour wave function of the nucleon for 
zero angular momentum is totally-symmetric, see App. 0. We identify 



TZ{p) 

nip) 









o r 




X 2 


® X 2 


® X 2 








i r 




x 2 ' 


®x^ 


® X 2 



*W(p) . 



18) 



One then obtains for the current matrix elements of the proton and neutron, respec- 
tively (see App. |y): 



(*p \iA %) 

( ^ n | J/i | ) 



1 

3 

+ 

1 

3 



#[°1 

^[1] 



■ diquark 
J u 

jdiquark 
J u 

■diquark 
J/j, 

• diquark 
J a 



+ 



3 \ 



■ quark 
J a 



^[0] \ 



• quark 
J a 



■ quark 
J u 



(19) 



Note that in the case of the proton, the photon coupling to the quark in the v-diquark 
channel drops out, whereas in the case of the neutron the current is sensitive to the 
difference between the quark and the diquark coupling. 



IV. RESULTS AND DISCUSSION 

The values of the parameters used in our model are given in Tab.|. The quark 
and diquark masses are the same as in our previous work [|lj]. The mass of the v- 
diquark is chosen to be equal to the scalar diquark mass. Also, the quark-diquark 
coupling constants are chosen to be equal: g = g s = g v . Their value is fixed by the 
minimum of the M(a) curve in the Ritz variational principle, analogous to Q. For 
six radial basis functions we have a minimum at M{pt) = 939 MeV for the oscillator 
parameter a = 1.35 fm. With equal masses and couplings we are closest to the 
SU (2) sp i n ® SU (2) isospin limit which is only broken by the different types of coupling 
(~ 1 for scalar and ~ 7 M for the v-diquark) (Eqs. (|7|) and @). Indeed, we find 
in this most symmetric case a best description of all electromagnetic form factors. 
The parameter influencing the shape of the form factors mostly is the diquark form 
factor parameter A, see It is put equal for both diquark types and is chosen to 
give a best description of the electric form factor of the proton. We obtain A = 0.24 
fm, in agreement with the commonly used diquark size ||. The proton electric form 
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factor is shown in Fig. ^ for momentum transfers 

up to -q 2 = 3 GeV 2 . We find 
a very good agreement with the experimentally found dipole shape. However, for 
about — q 2 > 3.5 GeV 2 , G P E becomes negative, though with small absolute value and 
converging to zero. This change of sign is the case for all four form factors and is 
model inherent. From the slope at q 2 = we get for the rms-radius of the proton 



(r 2 ) p = 0.84 fm, (20) 



compared to the experimental \J(r 2 ) p = (0.862±0.012) fm [19]. In Fig. |5| the neutron 

electric form factor is shown for momentum transfers up to —q 2 = 0.75 GeV 2 where 
is well established experimentally. We obtain a remarkable description of the 
experimental data. For the mean square charge radius of the neutron we find 

( r 2 ) n = -0.118 fm 2 , (21) 



which fits the experimental (r 2 ) n = -0.119±0.002 fm 2 |E0 . Since the neutron current 



is the difference of currents (see Eq. (0)), it is very sensitive to the parameters. The 
dashed curve shows G\ for a v-diquark mass m v = 670 MeV. Indeed, we describe 
the data well in the mass-symmetric case. In Fig. || we compare our results (full 
curve) with our recent calculation only involving scalar diquarks IJ (dotted curve). 
The improvement is obvious. 

Fig. [f] shows the magnetic form factor of the proton. As can be seen from Eq. 
( |17c| ), the spin flip current of the v-diquark is proportional to (1 + k), k being the 
anomalous magnetic moment of the v-diquark. This is chosen to be k = 1.6 to fit the 
experimental data for G P M . Again, our calculation reproduces the dipole shape very 
nicely. The dependence of G P M on k is shown by the dashed curves. The upper dashed 
curve corresponds to k — 2.0, the lower one to k — 1.0. The larger the diquark's 
anomalous magnetic moment the larger is the absolute value of the magnetic form 
factor. The neutron magnetic form factor is shown in Fig. ||. The calculation matches 
the experimental data qualitatively. However, the absolute value of the calculated 
curve is too small. Extrapolating to q 2 = 0, we get from G^(0) for the magnetic 
moments of the proton and neutron: 



(22) 



The SU(6) quark model prediction is — |. 

In Fig. [|we show the zero-components of the various currents (see Fig. [|). As shown 
analytically in |I[ , the corresponding quark and diquark currents have to be equal at 
q 2 = 0, thus equally contributing to the normalization of the Salpeter amplitude (Eq. 
(||)). It is interesting to compare these contributions of the three channels (scalar : 
zero-component of v-diquark : vector-component of v-diquark): 



Hp = 


2.78 /i at (exp.: 


2.793 fj, N ) 


Hn 


— 1.51 (exp.: 


- 1.913 n N ) 




-0.543 (exp.: 


- 0.685) . 
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0.63 : -0.23 : 0.60 , 



(23) 



which reproduces the result of |T2[ in the mass-symmetric case (where both compo- 
nents of the v-diquark-channel are added): 0.64 : 0.36. 

The four spin flip currents are shown in Fig. 10. It is mainly the current 



diquark 
3+ 



j (lowest dashed curve) which improves our results for the mag- 
netic form factors compared to our recent work It is interesting to see that at 
— q 2 > 1 GeV 2 the spin flip current, i.e. the magnetic form factor, is essentially due 
to the spin flip of the v-diquark. The transitions between the zero-component and 
the vector-component of the v-diquark, corresponding to the non-diagonal elements 
in the coupling matrices of the Eqs. ( |17a| )- (|17cj ) are seen to add up to zero: 



• diquark 
J u 



[1] 



[1] 



■diquark 
J a 



(24) 



see the dotted curves in Fig. |9] for /i = 0. 



V. SUMMARY AND OUTLOOK 

We developed a formally covariant constituent quark-diquark model of the nu- 
cleon. Two kinds of diquarks, scalar and axial-vector ones, are taken into account. 
The only interaction considered is a quark exchange. Starting from the Bethe- 
Salpeter equation for a bound state of a quark and a massive boson, we otain in the 
instantaneous approximation for the interaction kernel a pair of coupled Salpeter- 
equations. As a further approximation we choose a scalar-vector symmetric form 
of the spin-1 propagator. The equations are solved involving the Ritz variational 
principle with a finite number of radial basis functions. To obtain a stable solution, 
a form factor of the diquark has to be introduced. The current matrix elements are 
calculated in the Mandelstam formalism. With five parameters in the scalar-vector 
symmetric case (constituent quark and diquark masses, diquark size parameter A, 
diquark-quark coupling g and the anomalous magnetic moment of the axial-vector 
diquark) we find an excellent agreement with the experimental electromagnetic form 
factors for momentum transfers up to —3 (GeV/c) 2 . Only the magnetic form factor 
of the neutron is a little off the experimental curve. The calculated charge radii agree 
with the experimental ones, and we find the correct anomalous magnetic moment of 
the proton. 

It is indeed surprising to find such a good correspondance with the experimental 
data in such a simple model. To what extent this is due to the chosen quark-exchange 
interaction, which is essentially ~ - (~ ^ in coordinate space), remains to be exam- 
ined. One important ingredient is the covariant treatment and the correct boosting 
of the outgoing amplitudes. In what sense our results may provide a revealing of a 
probable diquark-structure in the nucleon is unclear. In any case, by introducing an 
effective potential and diquark parameters, the complicated three-body problem can 
be reduced to an effective two-body one. 
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A further important test of our model is the calculation of N-A transition form 
factors. In our model only the axial-vector diquark channel will contribute to this 
process. So the involved parameters should be further constrained. 

Acknowledgements: I am grateful to H.R. Petry, B.C. Metsch and C.R. Miinz 
for many helpful discussions and T. Schirmer for reading the manuscript. This work 
was supported by the Graduiertenkolleg 'Die Erforschung subnuklearer Strukturen 
der Materie'. 
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TABLES 



m q 


m s 


m v 


9 


A 


K 


350 MeV/c 2 


650 MeV/c 2 


650 MeV/c 2 


8.10 


0.24 fm 


1.6 



TABLE I. The parameters of the model 
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FIGURES 




FIG. 1. The (Bethe-)Salpeter equations for a bound state of a quark and a scalar or 
axial- vector diquark with a quark exchange interaction, see Eq. The dashed line 

represents a scalar diquark and the double-lined an axial-vector diquark. 




FIG. 2. The electromagnetic current is the sum of the diquark currents and the quark 
currents. 




FIG. 3. The coupling of the photon to a massive vector particle (Eq. ([Oj)). 
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G 



E 0.05 - 




0.4 

-q 2 [GeV 2 ] 



FIG. 5. The electric form factor of the neutron G^(q 2 ); experimental data are taken 
The full line corresponds to the parameter set of Tab. |, the dashed line 
670 MeV. 



from [19,21 
corresponds to m v 



• Platchkov et al. 
» Simon et al. 
■* Hanson et al. 
► Bartel et al. 




-q 2 [GeV 21 



FIG. 6. The electric form factor of the neutron G^(q 2 ); experimental data are taken 
from [19,23-25]. The full line corresponds to the parameter set of Tab. Q, the dotted line 
is our calculation with scalar diquarks only (from O). 
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1-0 , , 2.0 

-q [GeV ] 



3.0 



FIG. 7. The magnetic form factor of the proton G P M (q 2 ); experimental data are taken 
from [21,22]. The full line corresponds to the parameter set of Tab. fli the upper and lower 
dashed curve to k = 2.0 and k = 1.0, respectively. 



G 



-0.5 - 



M -1.0 




-2.0 



0.0 



• Hanson et a\. 

* Lung et a\. 

•> Hughes et a\. 
» Bruins et al. 



1-0 , , 2.0 

-q [GeV ] 



3.0 



FIG. 8. The magnetic form factor G 1 l I (q 2 ) of the neutron; experimental data are taken 
from |4||(]-!j§. 




1.0 

-q 2 [GeV 2 ] 
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FIG. 9. The charge density jo(Q 2 ) of the four currents of Fig. || in dimensionless units. 
The full lines describe the coupling to the quark, the dashed lines that to the diquark. 
The upper pair of full and dashed lines corresponds to the currents in the scalar diquark 
channel, the lower one to those in the v-diquark channel (with zero- and vector components 
added) . The dotted lines correspond to the non-diagonal couplings of Eq. ( |17a| ) . 



1.0 





\ 

- \ 

\ 

\ 

~x 




quark 

diquark 



0.0 1.0 2.0 

-q [GeV ] 

FIG. 10. The spin-flip currents j + = + i 22) (see Fig. ||) in dimensionless units. 
The full lines describe the coupling to the quark, the dashed lines that to the diquark. The 
full line starting with negative values is the quark current in the scalar diquark channel. 
The j+ current of the scalar diquark is trivially zero. 



APPENDIX A: PARITY TRANSFORMATION 



An axial- vector field transforms under parity-operation like [29 



VA»{x)V = -A^x) 



x M = x. 



Then, we find for the transformation of the Bethe-Salpeter amplitude: 

Xp(p)v = -Kv 7° Xp(PT ■ 



(Al) 



(A2) 



APPENDIX B: LORENTZ TRANSFORMATION 

For the Lorentz transformation of the v-diquark-quark Bethe-Salpeter amplitude 



we have 29 



However, as a spinor field, a nucleon field transforms like 
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*p(p) = Sa 1 ^ap(Ap) • (B2) 
This correct transformation is achieved by the covariant coupling 

*p(p) := 7 M XpW^ • (B3) 

APPENDIX C: NUCLEON WAVE FUNCTION 



Following Eq. fll8|) we write for the proton and neutron wavefunctions (without 
the totally anti-symmetric colour wavefunction) 

%(p)s = M (*i°J(p) xlf A + XMs) 

®n(p)s = Ar(y [ Hp)xl% + ^ ] (p)xlt) , (ci) 

where the flavour functions are mixed (anti-) symmetric: 

^.uud [0] 

Xm a — AO 

XT S = -^(xf 1 - V2 X l\) , (C2) 



with e.g. Xo = -^(ud + du)u. 
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